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Thermodynamical features of Xe system is investigated as functions of temperature and freeze-out 
density in the frame of lattice gas model. The calculation shows different temperature dependence 
' of physical observables at different freeze-out density. In this case, the critical temperature when the 

phase transition takes place depends on the freeze-out density. However, a unique critical excitation 
energy reveals regardless of freeze-out density when the excitation energy is used as a variable 
insteading of temperature. Moreover, the different behavior of other physical observables with 
, temperature due to different pj vanishes when excitation energy replaces temperature. It indicates 

' that the excitation energy can be seen as a more basic quantity to control nuclear disassembly. 

(N 

■ The phase transition and critical phenomenon of small systems is an interesting subject in recent nuclear physics 
" research. The break-up of nuclei due to violent collisions into several intermediate mass fragments (IMF), can be 
. viewed as critical phenomenon as observed in fluid, atomic, and other systems. It prompts the possible signature 
on the liquid gas phase transition of the nuclear system. On one hand, the onset of the multifragmentation [Q and 
vaporization ||^ channels can be seen as the signature of the boundaries of phase mixture . This is supported further 
by the fact that the caloric curve in a certain excitation energy range Q shows a saturate similar to a first order phase 
transition, in the framework of statistical equilibrium models On other hand, the observation of critical exponents 
parameters in the charged or mass distribution of the multifragmentation system Q can be interpreted as an evidence 
of the phase transition. Recently, the lattice gas model (LGM) has been applied to treat phase transition and critical 
phenomenon in the nuclear disassembly for isospin symmetrical and asymmetrical nuclear systems. LGM 
^ ' assumes a freeze-out density pf with thermal equilibrium at temperature T. The temperature was adopted naturally 
(~| ■ as a variable to study the feature of disassembly in nearly all previous calculation of LGM. In this paper, we will 
illustrate that the excitation energy can be taken as a more basic quantity to control the disassembly of nuclear system 
rather than temperature via studying the features of critical phenomenon and other physical observables in the lattice 
gas model. 

?H ' In the lattice gas model, A nucleons with an occupation number s which is defined as s = 1 (-1) for a proton 
(neutron) or s = for a vacancy, are placed in the L sites of lattice. Nucleons in the nearest neighbouring sites 

have interaction with an energy eg^sj- The hamiltonian is written hy E — 2m ~ '^i<j ^siSjSiSj. The interaction 

constant es^sj is related to the binding energy of the nuclei. Here enn,pp = = MeV, epn,np = = - 

5.33 MeV is used. The freeze-out density of disassembling system is pf = j^po where po is the normal nucleon density. 
The disassembly of the system is to be calulated at p/, beyond which nucleons are too far apart to interact. N + Z 
nucleons are put in L cubes with size by Monte Carlo sampling using the Metropolis algorithm. Once the nucleons 
have been placed, their momentum is generated by a Monte Carlo sampling of Maxwell Boltzmann distribution. 
Various observables can be calculated in a straightforward fashion. 

One of the basic measurable quantities is the distribution of fragment mass. In this LGM, two neighboring nucleons 
are viewed to be in the same fragment if their relative kinetic energy is insufficient to overcome the attractive bond: 
P^/2/i -|- e„p < 0. Once the fragment mass distribution is built, we can extract the effective power law parameters via 
fitting the mass distribution of fragments with Y{Ai) oc A^'^ and its second moment of fragment distribution defined 



X 



* Corresponding Author. Email: mayugang@publicl.sta.nct.cn 
^Mailing address 



1 



as JlOl 5*2 = — ^itAmai^ ^ where nj(Ai) is the number of fragments with Ai nucleons and the sum excludes the 

largest cluster A^ax- There are a minimum of r and a maximum of S2 at critical point for an infinite system. Besides 
the above quantities, we will use the average multiplicity < IM F > oi IMF and the information entropy H to search 
the critical point ||ri|,|l^ . H was defined firstly by Shannon in information theory ]l3t and can be introduced into 
nuclear dissociation ll^ , it reads H = ~ '^^Pi * ln{pi), where pi is the probability having "i" produced particles in 
each event, the sum is taken over all multiplicities of products from the disassembling system. H reflects the capacity 
of the information or the extent of disorder. 

We choose the medium size nuclei ^^^Xe as an example to analyze the nuclear disassembly. Three freeze-out 
densities of O.lSpo, 0.38poj and 0.60po, corresponding to the lattice size of 9^, 7^ and 6'^ respectively , were used. 
The calculations were performed from 3 MeV to 7 MeV and 1000 events were accumulated at each temperature and 
freeze-out density. 

We show the temperature and freeze-out density dependences of t, < IMF >, H and 5*2 in the left column 
of Fig.l. First, the critical temperatures determined by the extreme values are the same for the same freeze-out 
density, indicating that critical phenomenon really exsits. Second, the freeze-out density determines the temperature 
where the critical point is reached in this model. The critical point takes place at 4.25 MeV for freeze-out density 
0.18 po, 5.5 MeV for 0.38 po, and 6.5 MeV for 0.60 po, respectively. Obviously, the higher the freeze-out density, 
the higher the critical temperature. In this case, the critical behavior is determined by the two variables, namely 
temperature T and freeze-out density p/ as observed in Pan and Das Gupta [Q. But the extraction of T and pf 
is difficulty in view of experiments. On one side, the different thermometers give different apparant temperatures 
0,^. On other side, the freeze-out density is not directly measurable quantity in experiments. Moreover, the 
different initial conditions of models maybe result in different freeze-out density during the confrontations of model 
calculations with experimental data. All these facts complicate the accurate determination on temperature and freeze- 
out density and interfer with the correct extraction of physics. So it is interesting to search other variables to locate 
the critical phenomenon. A natural ideal is to adopt the excitation energy per nucleon E* / A, which can be deduced 
or reconstructed from the experiments, especially in At: detector measurements [ p^ . In the lattice gas model it can 

be defined as E* /A — Et ~ ^g.s. = + en,p—x^ — ^n,p—^, where N^p and N^]p is the number of the bonds of 
unlike nucleons at T and in the ground state, respectively. In a pure classical model the ground state corresponds 
to a cold nucleus at zero temperature and normal nuclear density where there is no kinetic energy and so that the 

ground state energy per nucleon is —en^p—^. Practically, N-(-_p is determined by the geometry and is equal to the 
maximum bond number of unlike nucleons possible for ^^^Xe. The Fig.le shows the excitation energy per nucleon 
E* /A in different temperature and freeze-out density. Noted that the curves of excitation energy are not linear with 
temperature. If performing the differentials for these curves, we can obtain an important thermodymical quantity: 
the specific heat per nucleon at constant volume (or density) as Ctj/A = The Fig. If shows C^/A as a function 

of T for the disassembling system ^^'^Xe at three p/. Clearly, the peaks of C^/A exist for the systems at each pf and 
locate closely at theirself critical temperatures, which supports strongly the viewpoint about critical feature as said 
above. By the mapping from E* /A to T, we replot the r, < IMF >, i7, and 52 as a function of excitation energy 
instead of temperature in the right column of Fig.l. Now nearly all the critical points locate at the same excitation 
energy regardless of the freeze-out density. Hence the excitation energy can be viewed as a more basic quantity in 
controlling the reaction dissociation. 

In order to illustrate this point further. Fig. 2 gives the average mass of the largest fragment {Amax): the isotopic 
ratio R(p/d) between protons and deutrons, the isobaric ratio R(t/'^iJe) between tritons and helium-3 and the ratio 
R(n/p) of emitted neutrons to protons as a function of temperature or excitation energy, respectively, in different 
freeze-out density. Again, the discrepancies stemming from different freeze-out density minimize when the excitation 
energy is used as the variable, and the curves for different pf merge approximately into a single line. 

To summarize, we studied the critical feature of Xe system and found that the critical temperature where the phase 
transition occurs changes with the freeze-out density which complicates the confrontation the theoretical predication 
with the experimental results. Contrary, we found that a unique critical excitation energy and the same excitation 
energy dependence of other physical observables reveals regardless of the freeze-out density, which indicates that 
the excitation energy can be viewed as a basic parameter to control the dissociation of the nuclear system. Unlike 
the temperature and freeze-out density, the excitation energy can be well extracted from experiments, especially for 
experiments using 47r multidetectors nowadays. Hence the use of excitation energy as a basic parameter will make it 
easier and definite to extract physics from the direct comparison between the experimental data and the theoretical 
calculation. 
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Figure Captions 

Fig.l: The observables as a function of temperature (the left columu) or excitation energy (the right columu) in 
different freeze-out density: the r parameter from the power law fit to mass distribution (a,g), the average multiplicity 
of intermediate mass fragments < IMF > (b,h), the information entropy H (c,i) and the second moment S2 (d,j). 
The mapping from temperature to excitation energy is plotted in Fig.le and the specific heat is shown in Fig. If. 

Fig.2: The average mass of the largest fragment in each event Amax the isotopic ratio R(p/d) between the 

protons and deutrons (b,f), the isobaric ratio R(t/'^_ffe) between the tritons and the Helium-3 (c,g), and the ratio 
R(n/p) of neutrons and protons (d,h). The left columu is ploted versus the temperature T and the right columu 
versus the excitation energy. 
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Fig. 1 




Fig.2 



